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Abstract 

Despite its importance for modeling the homogeneous hot early universe very little 
is experimentally known about the magnitude of the reheating temperature, leaving an 
uncertainty of remarkable 18 orders of magnitude. In this paper we consider a general 
class of polynomial inflaton potentials up to fourth order. Employing a Monte Carlo scan 
and imposing theoretical and experimental constraints we derive a robust lower limit on 
the energy scale at the end of inflation, > 3 x 10 15 GeV for sizable tensor modes, 
r > ICC 3 . If the reheating phase is perturbative and matter dominated, this translates 
into a lower bound on the reheating temperature, yielding T r h > 3 x 10 8 (7 x 10 2 ) GeV 
for gravitational inflaton decay through a generic dimension five (six) operator. 
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1 Introduction 

The primordial fluctuations of the cosmic microwave background (CMB) provide a unique 
window to energy scales far beyond the energies reached at colliders. With the recent Planck 
data [1], experiments have reached a remarkable sensitivity in decoding this information. The 
theoretical interpretation is necessarily prone to be model-dependent, with much progress in 
interpreting the recent data in terms of numerous specific inflation models, or classes of 
models, cf. Refs. m for a non-exhaustive list of examples. In this paper, we attempt to 
derive constraints on the relevant energy scales of the very early universe - in particular on the 
energy scale of the vacuum energy at the end of inflation and on the reheating temperature 
- in a more model-independent way, covering in particular all single-field inflationary models 
which can be approximated by a polynomial of degree four and a wide range of perturbative 
reheating scenarios. 
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The highest observable energy scale in the history of our universe is given by the poten¬ 
tial energy driving inflation at the time when the largest scales observable today exited the 
horizon, 

V* 1/4 = 1 r 1/4 M P ~ 3.3 x 10 16 GeVr 1/4 , (1) 

where A s denotes the amplitude of the scalar perturbations generated during inflation and 
Mp is the reduced Planck mass, Mp = 2.43 x 10 18 GeV. The unknown quantity is the tensor- 
to-scalar ratio r of the primordial fluctuations. Currently, the most stringent bounds on the 
tensor-to-scalar ratio r are obtained from a joint analysis of the BICEP2, Keck Array and 
Planck data, yielding r < 0.12 [4]. Many large-field inflation models predict values of r close 
to this bound, and upcoming experiments will probe values of r in the percent region, cf. 
e.g. [SHU- In this paper we focus on r > 10 3 , which might be probed in the near future. 

During inflation the potential energy decreases as the inflaton field slowly rolls down its 
scalar potential until finally the potential energy density becomes a subdominant component 
of the total energy density. This is the onset of the reheating phase, which produces a ther¬ 
mal bath with some high temperature T r p - the initial condition for the homogeneous hot 
early universe. The value of the reheating temperature is for example decisive for the ques¬ 
tion if thermal leptogenesis [8j can be responsible for the baryon asymmetry of our universe, 
which is only possible for T r h > 10 9 GeV [9]. Moreover, the reheating temperature controls 
the relic abundances of long-lived particles such as gravitinos present in supergravity theo¬ 
ries, which, depending on their mass, life-time and abundance, have the potential to explain 
dark matter [10, llj or to cause serious problems SMZ] for the cosmological evolution of 
the universe, in particular endangering the success of big bang nucleosynthesis (BBN). De¬ 
pending on the particle spectrum, a number of bounds have been derived on the reheating 
temperature mm- For example, in the context of the 17-parameter 72-parity conserv¬ 
ing phenomenological minimal supersymmetric standard model, gravitino dark matter with 
a stau as next-to-lightest supersymmetric particle is viable for a reheating temperature of 
T v h < 2 x 10 9 GeV |23j taking into account current collider bounds on the particle spectrum. 
However, despite its importance for modeling the very early universe, very little is experimen¬ 
tally known about the magnitude of the reheating temperature. It is bounded from below 
by the requirement of successful BBN and from above by the upper bound on the scale of 
inflation, 10 MeV < T r h < 2 x 10 16 GeV, leaving an uncertainty of remarkable 18 orders of 
magnitude. 

With the unprecedented sensitivity of CMB measurements, a number of authors have 
recently revisited constraints on the energy scales of the early universe. In a data-driven 
approach, the scalar potential can be reconstructed around the pivot-scale of a given experi¬ 
ment |24ll25j , cf. [Tj for a recent analysis. This yields rather model-independent information 
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on the scalar potential, however only in the observable range, i.e. for a few Hubble times or 
correspondingly for a sub-planckian range of the inflaton field. Constraints on the reheating 
temperature have been derived for given inflation models by extracting constraints on the 
number of e-folds elapsed after the pivot-scale left the horizon from CMB measurements [26]. 
Since the total number of e-folds elapsed depends (weakly) on the reheating temperature, 
matching the temperature and Hubble horizon today implies bounds on the reheating tem¬ 
perature which can be as stringent as T v h > 10 8 -10 10 GeV t23, though highly dependent on 
the underlying inflation model as well as on the equation of state during reheating [28j . 

In this paper, we pursue a different approach, aiming to reduce in particular the depen¬ 
dence on the choice of the inflation model. Performing a dedicated Monte Carlo scan based 
on polynomial inflaton potentials up to degree four and taking into account the recent Planck 
constraints on the inflationary observables, we obtain a lower bound on the energy density 
at the end of inflation in the single-field slow-roll approximation, > 3 x 10 15 GeV for 

r > 10 -3 - remarkably close to the scale of grand unified theories (GUTs). If the reheating 
phase is matter dominated, this translates into a lower bound on the reheating temperature, 
T r h > 3 x 10 8 (7 x 10 2 ) GeV for gravitational inflaton decay through a generic dimension five 
(six) operator. Both the bounds on I4 n d and T r h are strengthened for larger values of r. 

The remainder of this paper is organized as follows. In Sec.[2]we present analytical bounds 
on the energy density at the end of inflation and on the reheating temperature. These turn 
out to be highly dependent on the inflation model, motivating a systematic investigation of 
a broad class of models based on a Monte Carlo scan in Sec. [31 The results and discussion of 
this scan are presented in Sec. [I] before concluding in Sec. (5j 

2 Analytic estimates of the energy scales 14 n d and T r h 

Inflation is a well established concept to cure several short-comings of the standard model of 
cosmology, e.g. explaining the homogeneity and flatness of our universe on large scales [29] . 
Its predictions have been a striking success story, culminating in the measurement of the 
fluctuations of the CMB by the Planck satellite with remarkable precision, which strongly 
support the simplest model of single-field slow-roll inflation [I], Here, the dynamics of inflation 
is governed by a scalar field $ slowly moving down a scalar potential V (4>) , the latter providing 
a large, nearly constant energy density necessary to exponentially inflate the universe, 

3 = -V'itfi ), (2) 

with H denoting the Hubble parameter. Inflation ends at </> e nd = 4>(t e nt j) when the kinetic 
energy of the inflaton (f> becomes comparable to the vacuum energy V en ^ = V(0 en d)- Obser¬ 
vations of the CMB contain information about the scalar potential TV* Hubble times (e-folds) 
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earlier, when the pivot scale k* 1 2 of a given experiment exited the horizon, 

r ^end 

IV* = / “ Hdt ~ 50-60, (3) 

■It* 

corresponding to a field value of 0* = 0(t*) and an energy scale U* = V(</>*) 0 One of the main 
open questions in inflationary model building is the nature of the scalar potential between d* 
and </>end • 

By measuring the properties of the primordial fluctuations in the CMB, in particular the 
amplitude A s and the tensor-to-scalar ratio r, we can determine (or at least constrain) 14, 
cf. Eq. (H|). Can we use this information to constrain the energy scales at the end of inflation 
and reheating, 14 n d and T r h? After the pivot scale exited the horizon, inflation continued for 
a finite amount of time and the potential driving it cannot be arbitrary steep, so as not to 
violate slow-roll inflation. This will result in a lower bound on I4 n d> the value of the scalar 
potential at the end of inflation. This energy density is then converted into the energy density 
of the initial thermal bath. A crucial parameter is the time-scale of this transition, in the 
simplest setup governed by the decay rate of the massive inflaton particle after the end of 
inflation. Exploiting that the inflaton must at least interact gravitationally, the lower bound 
on 14 n d can be thus converted into a lower bound on the reheating temperature. In this 
section, our aim is to analytically examine the lower bound on I4nd( r ) which will then lead 
to a lower bound on the reheating temperature, T r h(r). 


2.1 Constraining the energy density at the end of inflation 


In this subsection we discuss analytical bounds on the drop of the potential energy during 
single-field slow-roll inflation. We discuss the most conservative case as well as the case of 
monomial potentials. These cases serve as reference cases for the later more general study 
presented in Secs. [3] and [4) 

Let us take two points in the inflationary period 1 and 2, t\ < t%. Using the potential 
slow-roll parameter^ 


e = 


mi (yi_ 
V 


Tj = Ml 


v ’ 


we consider the ratio of the potentials at point 1 and 2, 


Ui 

U 2 


= exp 


rVi 


'v 2 


dV 

~V~ 


= exp 


f r<t >i 


^ y/2e 
Mp 


d(j> 


(4) 


(5) 


1 More precisely, the CMB provides information on the shape of the scalar potential about 10 e-folds around 
this point. 

2 Effects of possible violation of the slow-roll approximation will be discussed when relevant. 
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In order to obtain a large value for V\jVi we require large e for as many e-folds, IV, as possible. 
On the other hand, in order to allow for large IV, e has to stay below e en( j (otherwise inflation 
ends). Hence, the largest V 1 /V 2 is obtained by a constant e just below e en( j for as many e-folds 
as allowed. For this case the number of e-folds is obtained by 


N = 


1 


r<t>i 

I(f) 2 \[2e Mp 


d (f> > 


4> 1 ~ <t> 2 


( 6 ) 


\/2e e nd Mp 

This determines the maximal field excursion allowed by slow-roll inflation. Hence, we can 
impose the conservative bound 

v, 

( 7 ) 


Yl p 2e end N 

V 2 


The corresponding potential yielding the largest 14/14 for a fixed number of e-folds has an 




exponential shape 3 


V(</>) = y ie ^7^i-^/ M P , 


( 8 ) 


and fulfills the condition p = 2e. Note that in this case, some other mechanism is needed to 
end slow-inflation (e.g. a second scalar held as in hybrid inflation). Requiring e, r/ < 1, Eq. ([7]) 
yields an absolute lower bound of 14/T4nd < e N *, corresponding to a decrease of five orders 
of magnitude in the energy scale for IV* = 50. However, as e ~ 1/2 is clearly excluded for 
the first 5-10 e-folds constrained by the CMB this bound is conservative and the dramatic 
exponential drop can at most take place for the last IV* e-folds not constrained by the CMB. 

Next, let us the consider the situation where 77 ( 4 ) = qe(4>), q E R. This corresponds to 
the case of a monomial potential V oc (jf, p > 0 m , and we can identify q = 2(p — 1 )/p. In 
this case 

V* = f £endV /2 = f 4e e ndA4 + P 

Fend \ e * J V P 

Note that Eq. Q contains the limiting case |ry| <C e for p ~ 1. Conversely, if \rj\ 2e we 
obtain 14/!4nd —>• 1, since the potential must be extremely hat. 

From this we can draw an important conclusion. In the special case of ry = 2e = constant, 
the amplitude of the scalar potential during the last IV* ~ 40 - 50 e-folds can indeed drop 
dramatically. However, this relies heavily on the specihc exponential shape of the potential 
and is not the ‘generic’ case. Indeed for monomial potentials, Eq. ([9]) yields a much more 
restrictive bound of (14/14nd) 1 ^ 3 4 < 5.2 for IV* < 60 and p < 3, whereby p > 3 is already 
strongly disfavored by the recent Planck data [IJ. In Sec. Owe will hence pursue the ansatz of a 
fourth order polynomial in order to systematically study the largest 14/14nd that additionally 
are in agreement with current observations from the CMB. In this ansatz it is impossible to 
achieve a constant value of rj = 2e for all values of (j) f° r non-vanishing ejf| In Sec. Owe will 

3 Models with an exponential inflaton potential were discussed e.g. in EDI- 

4 To maintain a constant value of q = 2e and hence saturate the bound Eq. J7| for IV* ~ 50 e-folds a 
polynomial of order ~ 180 is needed. 
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find that within this framework 14 /V erid is significantly smaller than the bound d7j) , namely 
14 /Vend < 10 -1000, depending on the value of the tensor-to-scalar ratio r. Moreover we will 
investigate the fine-tuning required in the parameters of the potential in order to saturate 
this bound. 

Before concluding this section, let us add a brief comment on the slow-roll approximation 
used in this section. To obtain large values of V*/V end we are interested in fairly steep 
potentials. In particular we will consider potentials which are well described by the slow-roll 
limit close to <£>* (as dictated by the Planck data) but may have sizable corrections as e ~ 1 
is approached towards the end of inflation. In other words, the violation of the slow-roll limit 
is typically accompanied by an accelerating inflaton field. Considering the full equation of 
motion for the inflaton, 

4> + 3Hj> = -V' ((/>), (10) 

we note that for monotonic potentials and an accelerating inflaton field, the slow-roll approx¬ 
imation over-estimates the velocity \(f>\ and hence the distance in field space \4>* — <^> e nd|- This 
effect is however typically (over)compensated since for an accelerating inflaton field, the slow- 
roll approximation under-estimates the value of |</>* — pencil when inflation ends, cf. footnote 0 
In general, the effect is however expected to be small since it only concerns the very end of 
inflation. We will return to this point in Sec. 14.21 

2.2 Linking inflation to the reheating phase 

Inflation is followed by a phase of reheating, during which the vacuum energy is converted 
into a hot thermal bath of standard model particles [32] . In the simplest setup, this implies 
an intermediate stage governed by the coherent oscillations of a heavy scalar field, which then 
decays into the standard model fields with a total decay rate T^ [33H35] , This field can but 
must not be the inflaton field itself. 

A scalar field oscillating in a potential V (0) = acfP implies an equation of state [36] 

= Oty) = P-2 
" <Po> P + 2’ 

where the energy density p$ = 4> 2 /2+V (< j>) and the pressure = <p 2 —V(cj>) are related through 
the virial theorem, p (V (</>)) = 2 (<fi 2 /2). Here we have assumed the oscillation frequency (j)/4> 
to be much larger than the expansion rate H, as is the case very soon after the end of inflation. 
Let us first assume for simplicity that the end of inflation is followed immediately by a matter 
dominated phase with uj = 0. In other words, we take the potential for the scalar held to be 
essentially quadratic: 

V(cf>) = ^m^((j) + c) 2 for 4> < 4> end with V ((/W) = Kmd > ( 12 ) 
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thus matching the energy densities at the end of inflation and at the onset of 
Finally imposing that inflation is indeed over at </> e nd@ 


reheating!^ 


eO) > e en d for 4> < (/’end 


(13) 


we find that the mass of the scalar field is bounded from below, 


ml > e end V end /Mp . 


(14) 


If indeed the inflaton field itself is responsible for reheating and the transition from inflation 
to reheating is perturbative, we can match not only the amplitude of the potential at </> e nd 
but also the first derivative. In this case the bound in Eq. (fT4|) becomes an equality. 

For completeness, let us consider also the case that the transition from inflation to a matter 
dominated phase is not instantaneous. For example, picture the more general situation that 
the potential for the scalar field </> is given by 


V($) 


—m 

2 


2 i 2 



w~ 2 

Mp ~ 2 



(15) 


for p > 2 where for simplicity we have chosen the origin of field-space so as V(0) = 0 = W(0). 
This corresponds to the situation where just after the end of inflation, while the inflaton 
oscillation is large, the potential is governed by some higher order polynomial until towards 
the end of reheating, as </> -C Mp, the quadratic term comes to dominate. According to 
Eq. (|11[) . this models an equation of state with uj in the interval {0,1}, in particular including 
the equation of state of radiation, u = 1/3, corresponding to a quartic potential. We do not 
consider 0 < p < 2 since the linear and constant terms in the scalar potential are always 
absorbed by our definition of the origin of field space. Fractional exponents in this range lead 
to potentials which are not differentiable at their minimum. 

Matching to V((/> en d) = Fend and e(</’ en d) > e e nd as above yields 



2 Fend 

6M| 



(16) 


5 Here, we are not specifying the details of the reheating mechanism. In particular, this may hold also for 
scenarios which involve tachyonic preheating m- In this case the large majority of the vacuum energy density 
is rapidly converted into low-momentum modes of a symmetry breaking field which again implies an equation 
of state with u> = 0 if the symmetry breaking field is sufficiently heavy. 

6 Note that more precisely, the end of inflation is reached when the first Hubble slow-roll parameter en = 
cj> 2 /(2H 2 ) reaches £h = 1 [38]. As long as slow-roll is a good approximation, ch — e. At the end of inflation, 
when slow-roll is violated, tH > 1, one finds for an accelerating inflaton field e(</>) > 1, allowing us the set 
tend = 1 in Eq. (1131) . In the (untypical) case of an decelerating inflaton field at the end of inflation, this only 
holds approximately. 
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where we have assumed that the (jp term dominates in the initial reheating stage (else it is 
irrelevant and Eq. (fTHl applies). We immediately see that for b ~ 1, i.e. if the quadratic and 
c jp term become comparable at <j) ~ Mp then Eq. (|16l) yields a comparable bound to Eq. (ED- 
In particular for p = 3 (p = 4) the bound on is weakened by merely a factor of about 2 (6). 
If b 2> 1, then the bound on disappears. This corresponds to the limiting case of a pure ()P 
potential. For b< 1, Eq. (11611 seems to improve on the bound of Eq. da, note however that 
in the derivation of Eq. (11611 we have assumed the <jp term to be dominant, i.e. b is bounded 
from below. We conclude that as long as the quadratic term is not strongly suppressed 
compared to higher order terms in the reheating potential, the bound in Eq. da is not 
significantly weakened. The following results are hence not very sensitive to the assumption 
of a purely quadratic reheating potential, or correspondingly to the assumption of u = 0. For 
concreteness, we will however focus on uj = 0 in the following. 

The reheating temperature is defined as the temperature of the thermal bath when half 
of the total energy density has been converted into radiation. This happens when the decay 
rate of the massive (^-particle becomes comparable to the Hubble rate, H ~ T^. Exploiting 
the expression for the energy density p r of a thermal bath at temperature T as well as the 
Friedmann equation, 


7r 


Pr = - g*T 

1 30 y 


Ptot — 3H 2 Mp , 


(17) 


where g* denotes the effective number of massless degrees of freedom in the thermal bath 


(<?* h = 427/4 for the standard model at high energies), this 


yield 


T r h = 


45 

7T 2 f£ h 


1/4 


yryiT] 


p ■ 


(18) 


A lower bound on this temperature can be obtained by assuming only gravitational interac¬ 
tions for the inflaton, resulting in the following effective dimension five operators: 

(19) 


\^4>b*b , , 

Mp Mp 


where B and / are bosonic and fermionic fields, respectively, and A is a dimensionless coupling. 
Such operators may be sourced by a coupling of the inflaton held to the kinetic terms of the 
respective particlesfl The decay rate of the non-relativistic particle cj) into a pair of light 


7 Hcre, we use H = F $ and p r = ptot/2 at T = T r h. Corrections arise due to the relativistic time delay 
in the decay of the ^-particle and from the precise factor a in p r = aptot|ir=r^,, which has to be determined 
numerically. However, even taking into account these effects, Eq. o typically remains a good estimate for 
the reheating temperature, cf. e.g. Ref. |39] . 

8 In the context of supergravity, such effective operators for the decay into a supermultiplet X are gen¬ 
erated e.g. by a term XMp 1 (/)\X\ 2 in the Kahler potential, which leads e.g. to a term XMp 1 c/>d l j,X*d fI 'X for 
the corresponding scalar field in the effective Lagrangian, see also m- However, determining the dominant 
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bosons or fermions is given by 

A 2 m? 

^ ~ 16vrM 2 ' 

In the case of multiple decay modes Eq. (1201) turns into the sum over the respective partial 
widths. One might also imagine the situation in which dimension-5 operators are forbidden 
or strongly suppressed. In this case there are additional suppression factors of m^/Mp with 
respect to Eq. (fl9l) . In summary, this yields 



T >dim5 

rh 


ndim6 

[ rh 


45 A 1/4 / \ 2 m z 6 


n 2 gl h 


16nMp ’ 


45 A 1/4 / A 2 m 3 


( 21 ) 


^r 2 gl h 


167rMp 


Plugging Eqs. (fT |) . (j7 |) , (fTTl) and (12TT) . this yields (for the case of the dimension five operator) 


T r d h im5 > 3.7MeV A 




2.2 x 10~ 9 , 

for the exponential potential and correspondingly with Eq. ([9]) for the monomial potentials 


T dim5 > 6.0 x 10 9 GeVA 


2.2 x 10" 9 


/ M 3 / 4 /427/4A 1/4 /£end\ 3 / 4 /(r,AA,e end ) 
vo.i/ V gl h ) vi/ o.o9 


(23) 

where f(r,N*,e end ) = (r/(16e eil d)) 3r ^* eend ^ 32eend_2r ^ (for typical values /(0.1,50,1) ~ 0.09). 
Eq. (1251) leads to T r h > 10 9 GeV for 0.001 < r < 0.1 and A ~ 1 with a very mild dependence 
on the choice of e en d and N *. On the contrary, in Eq. (|22j) these parameters enter in the 
exponent, implying that for N* e en d > 20 the bound does not impose any restrictions beyond 
the requirement that thermalization of the universe has to take place well before BBN, T r h > 
lOMeV. The huge difference between Eq. (1221) and Eq. (1231) stresses the need for a more 
model-independent approach, which we will address in the next section. 

Finally, we mention that other mechanisms can prolong the reheating phase by delaying 
the thermalization of the universe due to kinematic effects (such as e.g. kinematic blocking). 
These effects can further weaken the bound on T r h but are highly model-dependent mm- 


3 Bounds for polynomial potentials 

The previous section served to obtain an analytical estimate for a lower bound on the reheating 
temperature. In a next step, we will back this up by performing a systematic Monte Carlo scan 

contribution to the decay of the inflaton depends upon the SUSY breaking scale and the details of the super- 
symmetric particle spectrum. 
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of suitable inflationary potentials. We will specify the class of scalar potentials in Sec. 13.11 
Secs. 13.21 and 13.31 are dedicated to the imposed constraints. The procedure of the scan will be 
detailed in Sec. l3.4l while in Sec. l3.5l we will discuss the properties of the inflationary potentials 
particularly interesting for exploring the lower bound on the reheating temperature. Our 
results are then presented in Sec. [4j 


3.1 Generic inflaton potentials for sizable tensor modes 


Our aim in this section is to systematically sample a broad range of inflation models. We 
restrict ourselves to single-field slow-roll inflation models, implying that the potential must 
allow for approximately 50-60 e-folds of inflation and must allow for inflation to end, i.e. 
e(0end) ~ I- This in particular excludes situations in which the inflaton field becomes trapped 
in a false vacuum which it can only leave by tunneling through a potential barrier or if it is 
assisted by the dynamics of a further field. Instead, we require the potential to be sufficiently 
‘well-behaved’, i.e. to support a small but strictly monotonous slope over a sufficiently large 
field range. In this spirit, we expand the potential around the Planck pivot scale </>*, where 
it is well constrained by the CMB data, but truncate the expansion at fourth order: 


V(A<j>) ~ V* 



Mp d n V 
bn! d4> n 


(24) 


with Acf) = <fr — <j>*. The effect of the truncation is of course independent of the choice of 0*, 
e.g. all potentials of the type 

4 

V{<j>) = Y, (25) 

71=0 

are trivially included in the ansatz (1241) , including the monomial potentials with integer expo¬ 
nent p < 4. The motivation for this truncation (apart from the obvious effect of practicability) 
is manifold. For small-held inflation models, Acf> <C Mp, the Taylor expansion of any V((f>) 
around 4>* quickly converges and Eq. (1241) is a good approximation to the full potential. For 
large-held inflation models, which we are more interested in for the purpose of this paper since 
they come with a high scale of inflation, higher-order operators with sizable coefficients will 
make a sufficiently long period of slow-roll inflation difficult, unless these coefficients are very 
carefully tuned to ensure a cancellation of these contributions. Viable large-held models of 
inflation hence typically forbid these terms, e.g. by taking the inflaton to be a pseudo Nambu- 
Goldstone boson as in natural inflation }43li44| or, in supersymmetric models, by employing a 
shift-symmetric Kahler potential which is only broken by a renormalizable superpotential }45j . 
Our ansatz (j24[) thus covers a large range of known inflation models (cf. Ref. [I] for a recent 
overview in light of recent experimental results): Monomial potentials with non-integer power 
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p are well approximated by a polynomial of degree four, the crucial quantity V*/V en( \ differs 
by at most a factor of about 0.6-1.5 between the pure monomial and its approximation. 
It moreover trivially contains quadratic and quartic hilltop inflation [36], cf. Eq. (1251) . and 
approximately contains natural inflation, considering cos(</>//) (with i?!> <C /) truncated at 
fourth order. In a broader sense, it covers all renormalizable models of inflation. Finally note 
that typical small-r inflation models such hybrid inflation models [37[I38l|2| (putting aside the 
requirement of reaching e(</> en d) ~ 1) and Starobinsky-type inflation models [50] (including 
the original R 2 model as well recent rediscoveries of this exponentially flat potential [f6TU55j [) 
typically feature a very small I4/I4 n d < 10 and are hence not of particular interest for the 
lower bound on !4 n d and T r ^. 

Based on Eq. (1241) , we will systematically search for suitable inflation potentials by means 
of a Monte Carlo scan, restricting the coefficients of Eq. (1241) by the CMB observations and 
constraining the number of e-folds of inflation IV* by comparing the CMB pivot scale with the 
present Hubble horizon. We take the inflaton to be rolling towards the origin in field-space, 
i.e. A (f> < 0, and consider the range starting from A(f> = 0 and ending with e(Ac/> enc j) ~ 1. 


3.2 Constraints from Planck around A (j) = 0 


We consider the CMB observables A s ,n s , a s , n s ,r which can be expressed in terms of the 
potential slow-roll parameters e, r) and 


2 4 H'H (3) 

e = M £^- 


a 3 = Mi 6 


wV( 4 ) 

p y3 


(26) 


to leading order by 


A. = 


V 


24 AMj [1+0(e)1 ’ 


n s = 1 - 6e + 2p + O (e 2 ) , 
a s = lbery - 24e 2 - 2f + O (e 3 ) , 

k s = — 192e 3 + 192e 2 ry - 32 erj 2 - + 2pf + 2cr 3 + O (e 4 ) , 

r = 16e + O (e 2 ) , 


(27) 


where the slow-roll parameters are understood to be evaluated at Acfi = 0. For the computa¬ 
tions in our scan we use, however, expressions for the CMB observables beyond leading order 
which are given in the Appendix. We checked that in our scan the slow-roll parameters 

e =^ c m f? = 2c 2 , y/^ 2 = V 6cic 3 , \Zct3 = yj 24cf c 4 (28) 

9 These are listed under the name spontaneously broken supersymmetric models in [l]. Even a variant of 
the minimal hybrid inflation model constructed to achieve a large tensor-to-scalar ratio within an intermediate 
field excursion A0 en d < Mp only yields V,/V en d < 1.05 fUJ i. 
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are indeed sufficiently small compared to one in order to ensure the validity of the expansion 
(see e.g. [38]). This is in particular the case for potentials that allow for a sufficiently long 
period of inflation IV* > 50 and values of the potential drop of interest to our work, V"* /V en ^ > 
10, where the absolute values of the above parameters are found to be well below 0.1. 

The Planck collaboration has provided stringent constraints on these observables, which 
we take into account by calculating the % 2 of each parameter point v = (In A s , n s , r, a s ) T , 


x (A s ,n s ,r,a s ) = (v-vyC \v-v). 

Here v contains the mean values of the lcr confidence intervals measured by Planck 


(29) 


0 


ln(10 10 A s ) = 3.104 ± 0.035 , n s = 0.9644T, 0.0049 , 
r = 0 ± 0.069 , = -0.0085 ± 0.0076 , k s = 0.025 ± 0.013 , 


(30) 


and C is the corresponding covariance matrix, C = cr 1 Xa. Here cr is a vector containing 
the lcr uncertainties of the observables in Eq. (13011 and the correlation matrix X encodes the 
correlations between the observables. Note that here we assume that Acf) = 0 corresponds 
to the Planck pivot scale, a constraint we will explicitly impose in the next section. We use 
the data from the Planck legacy arxiv for the base-model including a non-vanishing running 
and tensor-to-scalar ratio, including the high multipole joint TT, TE and EE constraints 
and the low multipole polarization constraints [56]. This corresponds to confidence intervals 
given in the last column of Tab. 4 in the recent Planck analysis [T]. By reproducing the 
marginalized joined 68% and 95% C.L. constraints in the n s -a s - and n s -r- plane in mm we 
checked that the Gaussian approximation provides a good description of the uncertainties. At 
the time of writing, the Planck collaboration had not released an analysis based on a model 
including a non-vanishing r, a s and n s . We thus base our analysis on the data provided for 
the (r % 0, a s / 0, n s = 0) model. A priori, this could be a serious limitation since fairly 
strong correlations between these parameters have been pointed out, e.g. the a s -K s correlation 
was found to be important in [58] . However, performing the Monte Carlo scan detailed in 
Sec. 13.41 we find that the requirement of sustaining inflation for a suitable period (roughly 50 - 
65 e-folds, see Sec. l3.3l for details) enforces —2 x 10 -3 < n s < 10 , independently of the much 

weaker bound in Eq. (1301) . cf. left panel of Fig. |T] This justifies the use of Eq. (I29lh i.e. taking 
into account the correlations between A s ,n s ,a s and r but varying n s independently. In the 
following figures, we show the resulting 95% confidence region, corresponding to x 2 < 9.72. 


10 Of course the physical meaning of r as the tensor-to-scalar ratio implies r > 0 and we will consider only 
parameter points with r > 0. 
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3.3 Constraints on IV„ 


In the previous section we have focused on how the Planck data constrains the scalar potential 
locally, i.e. around A <f> = 0. Furthermore, the scalar potential must fulfill an important global 
property, it must sustain a suitable amount of inflation between A cf> = 0 and A</> enc j. This 
is ensured by comparing the CMB pivot scale fe* with today’s Hubble horizon. With the 
expansion history of the universe from reheating until today described by the standard model 
of cosmology, this determines the number of e-folds IV* of inflation elapsed after the scale 
fcp 1 exited the horizon: IV* = ln(ao/a*) + ln(a en d/ao)- Tracking the expansion history of the 
universe epoch by epoch, this yields (261159] : 


j\Tcond ^ 


66.64 


— InfJU+ilnf^ 


\a 0 H 0 J + 4 


Mp 


+ 


1 -3a; 


, 1 In (V* 

4 WendJ 12(1 + 0;) 


In 


( Prh 

V ^end 


(31) 


where u ; denotes the equation of state during reheating, p r h is the energy density related to 
the reheating temperature and fc* = a*I7* = 0.002 Mpc -1 is the Planck pivot scale, ao and 
Hq are the scale factor and Hubble constant today, where we work in units with ao = 1. 
The constant first term on the right-hand side is determined by the scale factor and Hubble 
parameter at matter-radiation equality ( a eq ,H cq ) and the degrees of freedom in the thermal 
bath at a eq and at reheating, gl q and g* h : 


In 


O'eqHeq 

aqHo 



3/-/ ( 2 q 

Mp 


1 a rh 

-In —(+ = 66.64. 

12 gt q 


(32) 


H eq and a eq are determined by the matter density today, = 0.3156, the red-shift of matter- 

radiation equality z eq = 3395 and the Hubble parameter today, Ho = h 100 km s _ 1 Mpc^ 1 , 
with h = 0.6727 [60]. As above g* is set to the corresponding standard model values, g* h = 
427/4 and g eq = 43/11. 

In the Monte Carlo scan we will impose the condition that IV* should lie within the interval 


N c°nd _ 3 < W < Arr d + 3 , 


(33) 


where we use ui = 0 as argued in Sec. 12.21 The reheating temperature is computed from 
Eq. (I2T1) and the corresponding inflaton mass is obtained from I7 en d according to = 
+ndEend/ Mp. The latter corresponds to Eq. (fl4l) if a single field 4> is responsible for driving 
inflation and reheating. If this is not the case, larger values of and hence p v h are possible. 
However IV* only depends very weakly on this quantity, AIV* = (In Ap r h)/12. For simplicity 
we hence use the equality instead of the more general inequality in Eq. (1141) for evaluating 
IV*, and take the uncertainty to be covered by allowing IV* to be in a finite range around 
jy-ccmd, c £_ gq. ([ 33 ]) . On the contrary, the lower bounds on I4 n d and T r h (cf. Eq. (|2IT) 1 are 
based on the more general inequality in Eq. m, hence the results we present in Sec. [4] are 


14 








absolute lower bounds. If the operator of the inflaton decay is not specified we impose the 
conservative condition 


iV“ nd (dim 6 ) — 3 < TV* < N™ nd (dim 4) + 3, (34) 

with I\T ond (dim 6 ) referring to Eq. (13111 with p r h expressed by means of the dimension six case 
in Eq. m • Accordingly, iV“ nd (dim 4) corresponds to the dimension four (i.e. unsuppressed) 
case, i.e. removing m^/Mp in Eq. (fT9l) . In both cases, we set the coupling constant A to one. 

3.4 Monte Carlo scan 

In this subsection we briefly summarize the computational steps of the Monte Carlo scan over 
the parameter space spanned by the model (1241) . i.e., the parameters 14,ci,C 2 ,C 3 ,C 4 . These 
parameters can be expressed as analytic functions of the CMB observables A s ,n s ,r,a s , k s to 
leading order in the slow-roll parameters (using (12711 1 in order to obtain an estimate of the 
parameter region of interest, which we use in the course of generating scan points as described 
below. However, after the generation of a parameter point we only keep the parameters 
14, ci, C 2 , C 3 , C 4 and recompute the CMB observables using the full expressions shown in the 
Appendix. 

For each generated point of the Monte Carlo scan we perform the following steps. We 
first randomly choose A s , n s using a Gaussian probability distribution around their mean 
values (see Eq. (13U14 with various widths up to 3 standard deviations as well as r with 
a logarithmically flat distribution between 10 -3 and 0.3. This determines 14, ci and C 2 
(at leading order in the slow-roll parameters). For variation of the coefficients C 3 and C 4 , 
determined by a s and k s at leading order in the slow-roll parameters, we follow three different 
approaches used for generating different sets of parameter points. 

1. Randomly choosing a s and n s using a Gaussian probability distribution around their 
mean values (see Eq. m) with up to 2 standard deviations. 

2. Randomly choosing C 3 , C 4 using a logarithmically flat probability distribution with dif¬ 
ferent ranges between 10 _1 and 10 -15 . 

3. Randomly choosing C 3 , C 4 close to values that lead to a vanishing discriminant (and 
discriminant of the discriminant) of the potential (see Sec. 13.51 for a detailed discussion). 

The complete set of parameter points has been found to provide a sufficiently dense coverage 
of the region of large potential drops 14 /I4nd and therefore ensures the formation of a sharp 
transition between forbidden and allowed points. Note that the absolute density of scan 
points does not claim any physical meaning. 
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Figure 1: Scan points in terms of CMB observables. Left panel: Running (a s ) and running of the running 
(k s ) of the spectral index. Right panel: Spectral index n s and its running a 3 . Color-code: Blue: Inflationary 
points. Yellow: Inflationary points within 2-sigma-contours regarding Planck data at A <f> = 0. Green: Points 
additionally fulfilling the constraint on TV*. 



For each point we then compute A^ en d, requiring 


e(A(jiend) — ^end j 


(35) 


where we use e en( j = 1 (later we will also consider e enc ] = 3 for an estimation of the systematic 
uncertainties due to the use of the slow-roll approximation, cf. also footnote [6]) . The solution 
of Eq. (1351) only corresponds to a physical solution, if A^ en d < 0, V'(A4> en d) > 0 and if there 
is no zero of V' between 0 and e en d- Only points that fulfill these requirements are kept as 
potential inflationary points. For each potential inflationary point we compute IV* from 


A* = 


r° V d(A0) 


(36) 


Points with 10 < A* < 10 5 are kept. Further we compute the y 2 as described in Sec. 13.21 and 
impose the consistency condition for A* arising from the comparison of the CMB pivot scale 
with the present Hubble horizon as detailed in Sec. 13.31 We generate a total of 15 million 
inflationary points. 

Fig. Q] shows our scan points in the a s -K s and n s -a s plane. Inflationary points (i.e. poten¬ 
tials which support slow-roll inflation for 10 < A* < 10 5 e-folds) are marked in blue whilst 
points that additionally fulfill the Planck constraints on the scalar potential at A <j> = 0, as 
described in Sec. m are marked in yellow. Note that these points do not represent potentials 
fully compatible with the Planck data, as the Planck data additionally contains the informa¬ 
tion that Ac/) = 0 corresponds to a pivot scale of fc* = 0.002 Mpc -1 . Rather these yellow 
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points (referred to as “locally Planck allowed” in the following) represent potentials which 
are locally consistent with the Planck data around A<p = 0, but in general will not give the 
observed amount of inflation between A <f> = 0 and A0 enc j. This final condition is imposed by 
means of the consistency condition Eq. (1341) and results in the points marked in green. As 
mentioned in Sec. 13.21 points that fulfill Eq. (1341) only survive in a very small interval in k s . 
The stray blue points around a s ~ 0, n s < —0.02 in the left panel correspond to very small 
values of I4/T4 n d and low values of N .*, close to the lower bound N .* > 10. Since this region 
of the parameter space is irrelevant for the lower bounds on V en & and T r h, the density of scan 
points is very low here. 

3.5 Properties of the inflaton potential 

The parameter region of particular interest to this study is the regime which allows for the 
largest decrease of the potential energy during inflation. This implies a fairly steep potential 
allowing for a large energy drop during the final IV* e-folds while at the same time keeping 
the slow-roll parameters sufficiently small. In particular, not only must the derivatives of the 
potential be small compared to the inflationary scale 14 , but they must ‘track’ the rapidly 
falling energy density V (A (f>) until after A* e-folds the slow-roll parameter e reaches its critical 
value and inflation comes to an end. In Fig. [2] we schematically show the type of potentials 
which meet these requirements, as found in the Monte Carlo scan. Here, we have imposed 
the 95% confidence interval from the CMB data, cf. Sec. m and required values of A* 
in accordance with Eq. (1341) . The color-coding refers to the tensor-to-scalar ratio r. We 
immediately notice that we can identify two different regimes, corresponding to potentials 
with minimum (blue) or saddle point (green) close to the end of inflation. 

This can be understood as follows. As V (A cf>) decreases towards the end of inflation, the 
slow-roll parameter e increases and inflation comes to an end - unless the numerator V'(Acj)) 
of e vanishes simultaneously. The smallest values for I4nd are hence achieved if V(Acf>) and 
V'(A(f>) have a common zero point at a suitable value of A(j>. This implies an (at least) double 
zero point of V (A (f>) and hence a necessary condition for the smallest I 4 n d is a vanishing 
discriminant of V(Acp). In practice, small values of I4 n d after a finite amount of e-folds 
are achieved if this condition is approximately satisfied, i.e. if the zero-points approximately 
coincide, 

Disc[F(A0)] ~ 0 . (37) 

Eq. (1371) is an equation of degree three for the parameter C 4 with generically one real and two 
complex solutions. We observe that potentials with particularly small V en( i are achieved if (at 
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Figure 2: Scalar potentials which maximize the energy drop during the final IV* e-folds of inflation for different 
values of the tensor-to-scalar ratio r (color-coded on a logarithmic scale). The orange line marks the end of 
inflation at e = 1. 

least) two of these solutions approximately coincide, i.e. 

Disc[Disc[F(A^)](c 4 )] ~ 0 . (38) 

This condition implies that Disc[U(A0)] has a double zero point at C 4 € M, ensuring that 
the smallness of Disc[U(Ac)>)] is stable against small variations of the parameter C4, allowing 
us to find solutions for V(A(j>) with (several) close, but not exactly coinciding zero-points 
- precisely what we are after. The two distinct situations observed in Fig. [2] correspond to 
potentials with two (blue) or three (green) almost coinciding zero-points. Four coinciding 
zero-points would simply correspond to a </> 4 potential, which is however clearly excluded by 
the Planck data. Eqs. (1371) and (1381) provide us with analytical solutions for two parameters 
(04 and C3). Using values in the close vicinity of these solutions significantly increases the 
efficiency in exploring the large U*/ 14 n d region. 

Moreover, we can exploit that certain ranges of the coefficients c n lead to particularly 
large energy drops during inflation. For example, small values of the higher coefficients lead 
to well-behaved potentials, which are more likely to support inflation over a sufficient amount 
of e-folds. Also, we note that for small values of r, the A(ji 3 -coefficient is indirectly proportional 
to the square root of the tensor-to-scalar ratio, C3 oc —a. s /y/r oc 1/ci. This correlation is only 
strengthened by imposing the requirement on A*. Implementing these observations allows us 
to efficiently investigate the region of large U*/U e nd- 




r = 0.001 


r = 0.15 
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4 Results and discussion 


In this section, we present and discuss the results of the Monte Carlo scan. Following the 
strategy of the analytical discussion, cf. Sec. [2J we will first discuss the constraints on the 
scalar potential during the final N. * e-folds before drawing conclusions about the following 
reheating process. 

4.1 The energy density at the end of inflation 

The dynamics of inflation, described by Eq. ([3]), is independent of the overall scale 14 of 
the potential as defined in Eq. (041) . As a first step, let us hence consider the relative drop 
of the scalar potential during inflation, 14/Vend, as a function of N * and r, respectively, cf. 
Fig. [3j Here we mark in blue all parameter points which sustain inflation for IV* = 10-10 5 
e-folds - implying that these potentials also contain a mechanism to end inflation by reaching 
e(A0) = e end = 10 In yellow, we show the parameter points which additionally locally 
match the Planck data at 95% C.L., cf. Sec. 13.21 Finally in green, we impose the consistency 
condition on IV* given by Eq. (13411 . For comparison, the black dashed lines show the results 
obtained in Sec. [2] for the monomial potentials V oc cjp and the exponential potential V oc 
exp(v / 2e end A0/Mp), cf. Eqs. © and For monomials with p < 4, the red dashed lines 
indicate the range of fV*-values in accordance with Eq. (1341) . 

Focusing on the locally Planck allowed points in the left panel, we first note that all points 
covered by our ansatz and locally allowed by the Planck data are far below the absolute upper 
bound for 14/Vend given by the exponential potential. Moreover we can classify these points 
by the number of approximately coinciding zero-points at A <fi ~ A^ en d, cf. Sec. 13.51 Points 
with one, two or three approximately coinciding zero-points are predominantly located in 
three horizontal bands yielding increasing values for 14/14 nd , roughly bounded from below 
by 14/Eend = (1, 25, 250}, respectively, and broadening for increasing values of A*. In the 
lowest band E(A <j>) has a single (approximate) zero-point. As A <f> approaches this value, 
the slow-roll parameter e blows up rapidly and inflation ends. The higher bands correspond 
to two and three (approximately) coinciding zero-points, i.e. common zero points in V(Acp) 
and V'(A<f>), allowing e to keep a finite value as V(Acf>) approaches very small values. The 
monomial potentials with p = {1, 2, 3} can be seen as prototypical examples. 

Also in the r- 14/!4 n d plane, the points locally allowed by the Planck constraint can be 
grouped according to number of coinciding zero-points, yielding three horizontal bands which 
turn upwards at r ~ 0.1-0.2. The structure of the resulting yellow region in the right panel 

11 In this sense, we restrict our analysis to single-field inflation models, excluding in particular models of 
hybrid inflation where the dynamics of a second scalar field is responsible for ending inflation. 
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Figure 3: Scan points in the Ah-^/V/ndl-plane and r-(\ 4/14 n d)-plane. Color-code: Blue: Inflationary points. 
Yellow: Inflationary points within 2-sigma-contours regarding Planck data at A0 = 0. Green: Points ad¬ 
ditionally fulfilling the constraint on TV*. Dashed black curves: Prediction for monomial potentials 4> p with 
p = {1/2,1, 2, 3, 4} with the constraint on N * fulfilled within the region marked dashed red (for p < 4). In the 
left panel, we additionally show the most conservative analytical bound, the exponential potential. 


of Fig. [3] can be understood based on Eq. (|37l) . Taking a s and k s to be close to zero as 
indicated by the Planck data and n s within the 2a region, Disc[P(A</>), Ac/>] is approximately 
constant for small r until reaching a zero-point around r ~ 0.11-0.23. In this region of 
large r, very small values of Disc[P(A</>), A</>] are reached, corresponding to (several) nearly 
coinciding zero-points allowing for large V*/V en d- In this sense, the recent Planck results have 
now reached a crucial sensitivity, excluding the part of the parameter space where V*/V enc \ 
begins to shoot up. Correspondingly, if future measurements are able to constrain the tensor- 
to-scalar ratio to be below about r ~ 0.05, the upper bound on 14/V/ n d will become about 
an order of magnitude stricter. The overall shape of the green region can be understood by 
employing monomial potentials as prototypical examples. High values of 14/I4nd are reached 
if V(Acp) has a multiple zero-point (j> > 2, corresponding to a sizable value of r) whereas at 
small values of r, for p <C 1, the potential becomes very flat and I4/I4 n d —> 1- 

From the discussion above we can anticipate that the parameter points with high 14/Vend 
come with finely tuned parameters. To quantify this, we determine the fine-tuning mea¬ 
sure [6TJ 


£ = Max 


S(V,/V e nd) 


, n = 1,...,4 


(39) 


5c n (K/Pend) 

with 5c n denoting a small variation of the respective parameter of the scalar potential 
\5c n /c n \ ~ 10- 6 . Large values of £ indicate a high degree of fine-tuning, implying that the 
value of i4/{4 nd is extremely sensitive to at least one of the parameters. On the other hand, 
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Figure 4: Fine-tuning £ among the locally Planck allowed points in the r-(V*/V en d)-plane (left) and lower 
bound on the energy density at the end of inflation in the r-14 n d-plane (right). The color-code in the right 
panel is as in Fig. [3] 

points with £ = 0(1-10 ) correspond to potentials which do not require any particular tuning. 
In the left panel of Fig. [4] we show the fine-tuning measure £ where we performed a binning 
among the locally Planck allowed points in the r-(V* /h^nd) plane and displayed the point with 
the smallest value of £ in each bin. As expected, high values of V*/V en( \ always come with a 
high degree of tuning. This is rather independent of the value of r. In particular, the three 
horizontal ‘bands’ mentioned above, corresponding to an increasing number of coinciding 
zero-points, are also clearly distinguished by increasingly higher fine-tuning. Dropping the 
requirement of lying within the 2-sigma-region of the Planck data and considering the full 
data set, we observe that the fine-tuning is relaxed for small r by up to a factor of 10 but 
remains basically unchanged for r ~ 0.1. The high tuning is hence mainly a result of the 
careful balance of V and V' in the slow-roll parameter e, and is not a feature of the specific 
values of the CMB observables, in particular for large values of r. Comparing the left panel 
of Fig. [J] with the points passing all constraints in the right panel of Fig. [3l we note that in 
particular the large values of I4/I4nd around r ~ 0.1 feature a maybe uncomfortably high 
degree of tuning. A restriction to less severely tuned potential would correspondingly yield 
a tighter upper bound on lA/V^nd; translating to a tighter lower bound on the reheating 
temperature. Since choosing the acceptable degree of tuning is however a very subjective 
question, we continue in the following with the full data set. 

Finally, in the right panel of Fig. [U we include the information on 14 which comes with 
an additional r-dependence. This leads to the final constraints on the energy density at the 
end of inflation, one of the main results of this paper. Independent of any assumptions about 
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the reheating process which we will make in the following section, we find that the CMB 
observations constrain this energy scale to 3 x 10 15 GeV < (I4 n d) 1,/4 < 2 x 10 16 GeV for 
r > 10 -3 for the scalar potential of any single-field inflation model described by operators up 
to dimension four. This remarkable narrow range close the GUT scale might be taken as an 
indication to explore connections between inflation and the breaking of a unified gauge group. 
While this has been recently much discussed in the context of a large tensor-to-scalar ratio, 
we stress that indeed this conclusion holds over the entire range of the tensor-to-scalar ratios 
investigated, 0.001 < r < 0.15, since the overall scaling of the potential, V* oc r, is to some 
extent compensated by the large decrease in energy allowed by large r. This is illustrated by 
the example of the monomial potentials, cf. dashed red lines in Fig. |4j For small values of r, 
the drop in energy during inflation becomes negligible and the energy at the end of inflation 
is mainly determined by the r-dependence of V*. As an aside, we point out that this provides 
a good indication on how to extend our bounds to smaller values of r < 10 -3 . The results 
depicted in the right panel of Fig. [4] stress the importance of upcoming CMB experiments. 
The detection of a tensor-to-scalar ratio in the percent regime, 0.01 < r < 0.06, would tighten 
the bound on the energy density at the end of inflation to ( V end) 1 / 4 > 5 x 10 15 GeV. 

Bounds on the scalar potential of inflation have been derived in other contexts. Using the 
same ansatz as Eq. (1241) . the Planck collaboration recently performed a local reconstruction 
of the scalar potential around 0*, restricting themselves to the observable part of the poten¬ 
tial [TJ. We find that our ranges for the parameters c n lie well within the ranges found in this 
local reconstruction. This is an important cross-check, implying that also a few e-folds before 
the pivot-scale our potentials are sufficiently well behaved to not generate any anomalies at 
low multipoles in the CMB. Compared to Ref. [lj, our parameter ranges are of course much 
more constrained since we additionally impose the requirement that extending Eq. (|24j) to 
large values of A (j>, the scalar potential supports slow-roll inflation for a suitable amount of 
e-folds. Moreover, recently in Ref. [5Ej a lower bound on V*/!4 n d was derived, V*/V e nd > 1-01 
for r < 0.15, in this sense focussing on the opposite limit as the bound presented here. 

4.2 A lower bound on the reheating temperature 

Equipped with a lower bound on the energy scale at the end of inflation, I4nd, we can proceed 
to constrain the reheating temperature based on Eqs. (|14D and (1211) . In Fig. [5] we show the 
resulting lower bound on ^rh based on a dimension five and dimension six operator for the 
inflaton decay, respectively. Based on the monomials as prototypical examples the overall 
scale and shape of this constraint is well understood by applying Eq. (1231) . In both cases, we 
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obtain values significantly above the BBN constraint of about 10 MeV, 


T r d h im 5 > 3 x 10 8 GeV, 
T r d h im 6 > 7 x 10 2 GeV. 


( 40 ) 


for r > 0.001. Here, as in Fig. [5] we have set the coupling constant in Eq. (1211) to A = 1. Note 
that T r h depends linearly on A, hence a suppression of this coupling constant can weaken the 
bound accordingly. 




Figure 5: Scan points in the r-Tlh-plane showing the lower bound on the reheating temperature for a decay via 
dimension five operators (left) and dimension six operator (right) with A = 1. The color-code is as in Fig. [3] 

The left panel of Fig. [5] displays a remarkably strong bound for the dimension five case, 
implying that temperatures at least very close to the temperatures relevant for the gravitino 
problem and thermal leptogenesis were indeed reached in the early universe (T 10 9 GeV). 
This bound might be further tightened if upcoming B-rnode searches succeed in detecting a 
tensor-to-scalar ratio of a few percent, 0.01 < r < 0.07. In this region, the lower bound on 
the reheating temperature for A = 1 reaches as high as T r h > (1-2) x 10 9 GeV. This bound 
is however significantly weakened if the coupling A is suppressed, either by a symmetry or by 
an accidentally small value. To demonstrate this effect, we consider the situation where the 
dimension five operator is negligibly small and the time-scale of reheating is instead governed 
by a dimension six operator in the right panel of Fig. 0 In this case the constraint on the 
reheating temperature is significantly weaker, owing to the additional factor of m<p/Mp ~ 
y/V en d/Mp, but is still well above the BBN bound of 10 MeV. 

Throughout this study we have employed e = 1 as the condition for the end of inflation. 
Since slow-roll is violated at the end of inflation, this is however only approximately correct, 
more precisely inflation ends at en = cp 2 /(2H 2 ) = 1 [38]. The effect of this on the bounds 
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derived on the quantities of the reheating process is discussed in footnote El here we comment 
on the resulting uncertainty of the quantities V*/V eri( \ and T r h due to a misestimation of 
A0 enc p To this end we have solved the full second order differential equation m for a 
representative subset of 20000 parameter points^ (Unfortunately, applying this procedure 
in the full parameter scan is computationally too demanding.) As expected, we find e(en = 
1) > 1, hence the slow-roll approximation underestimates A(/> en d and thus U*/U e nd- We find 
that without the slow-roll approximation, A* increases by AA* = 1-1 - 2.2 and I 4 /I 4 n d 
increases by an 0(1) factor, which in the most interesting range 45 < A* < 65 is found to be 
(U*/Uend) full /(K/U en d)sr — 1-8 - 2.1. We observe that since e grows very rapidly at the end 
of inflation in the models of interest here, the errors in A^> en( j and the derived quantities are 
small. The bound (|14l) on m^, entering the bounds on the reheating temperature via Eq. (1211) . 
changes by a factor of 0.7 - 1.8. This moderate change is due to compensating effects of e en d, 
which increases, and !4 n d, which decreases with respect to the slow-roll approximation. In 
summary the uncertainty arising from the use of the slow-roll approximation is very moderate. 

In the analysis presented here, we have restricted ourselves to a matter-dominated reheat¬ 
ing phase, oj = 0, while allowing for a wide range of inflationary potentials^ An alternative 
approach to constrain the reheating temperature is to fix the inflation model and then use 
Eq. (1311) to constrain the reheating temperature. This allows one to probe all reheating sce¬ 
narios with a constant u. A recent analysis along these lines can be found in Ref. m- For 
simple inflationary models, such as the monomial potentials which have served as prototypical 
examples throughout this paper, the resulting bounds on the reheating temperature can be 
determined analytically. For a cj ) 2 monomial potential and cj = 0, the reheating temperature 
is found to be larger than about 10 9 GeV from the lcr constraints on the spectral index 
n s , consistent with our results. However, lowering the power of the monomial to p = 1, 
the bound completely disappears, allowing for reheating temperatures in the full range of 
10 MeV < T r h < 2 x 10 16 GeV. The result is also very sensitive to the choice of u, increasing 
its value to uj = 2/3 implies for the p = 2 case that no significant bounds can be placed on the 
reheating temperature, whereas for uj = — 1/3 the lower bound is raised to about 10 14 GeV. 
While it is remarkable that the CMB data alone can yield strong constraints on the reheating 
temperature for certain inflation models, these results also stress the need for a more model 
independent approach, as we have pursued in this paper. 

Further approaches to constrain the reheating temperature include constraints from the 

12 To obtain robust statements, we consider three different parameter sets: After binning our points in the 
r-(V*/l4nd) plane, we consider first a randomly chosen point within each bin and, second, the point with the 
highest and lowest fine-tuning, respectively. Third, we re-compute our results for e en d = 3 (which is found to 
be a typical value for e(tH = 1)), and select the allowed points with highest V*/V en d- 

13 Note that as argued around Eq. tm we expect the dependence of our analysis on the value ui to be small. 
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thermalization of neutrinos, measured by the effective number neutrinos iV e ff required to 
reproduce the CMB spectrum |62lf63l . implying however only a rather weak bound of T r h > 
3.2 MeV. Moreover, a detection of gravitational waves by future space-based experiments 
such as BBO [61] and DECIGO [65] would open a completely new window to the energy scales 
of the early universe. Given a sufficiently large tensor-to-scalar ratio, this would allow to probe 
the reheating in a range of 10 6 - 10 9 GeV [ 66 ]. Finally, from a theoretical point of view it has 
been argued that there is an upper bound on the reheating temperature, T < 10 9 GeV, since 
large couplings of the inflaton field to any other particles would destroy the flatness of the 
inflaton potential via radiative corrections [32]. This bound can however easily be avoided if 
the time-scale of reheating is set by the decay of a particle different than the inflaton, as can 
occur e.g. in hybrid inflation or in some preheating scenarios. 

5 Conclusion 

In this paper we derived constraints on the relevant energy scales of the very early universe, 
namely the energy scale at the end of inflation, I 4 n d, the mass of the inflaton, m and 
the reheating temperature, T r h. The most dramatic potential drop V*/V en & is realized in an 
exponential potential yielding the most conservative bound V*/I4nd < e N * ■ However, for a 
realistic single-field inflation scenario - providing a mechanism to end inflation as well as being 
consistent with observations of the CMB around the pivot scale - a much stronger bound can 
be achieved. In order to cover a large class of relevant inflation models we examined a general 
fourth order polynomial potential. We constrained this potential by requiring consistency 
with recent data from the Planck collaboration and by demanding that inflation ends after 
the number of e-folds which provides a matching between the CMB pivot scale and today’s 
Hubble horizon. We found that this restricts the potential drop to stay below I4/I4nd — 1000, 
yielding the highest values for large tensor-to-scalar-ratios close to being excluded by Planck, 
r ~ 0.1. For smaller r the potential drop becomes significantly smaller. This results in the 
bound > 3 x 10 15 GeV for r > 0.001 - remarkably close to the GUT scale. 

Matching the inflaton potential at the end of inflation to the subsequent oscillation phase 
allowed us to connect I4 n( j and which provides a lower bound on the inflaton mass. This 
allows us in turn to estimate the decay width of the inflaton governing the time scale of the 
energy transition to the thermal bath and hence the energy density at the beginning of the 
radiation dominated phase. Without taking into account the decay width of the inflaton, 
the reheating temperature can only be constrained by solving Eq. (1311) for p^. However, due 
to the logarithmic dependence this procedure is subject to huge uncertainties, in particular 
it does not provide an interesting bound in the general class of models considered in this 
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paper. Instead we specified the decay mode of the inflaton by exemplarily considering generic 
Planck-suppressed dimension five and six operators. This allowed us to place bounds on 
the reheating temperature, > 3 x 10 8 GeV and T^ 1 ™ 6 > 7 x 10 2 GeV, significantly 

above the requirement that thermalization of the universe has to take place well before BBN, 
T r h > lOMeV. If moreover upcoming B-rnode searches succeed in detecting a tensor-to- 
scalar ratio of a few percent, 0.01 < r < 0.07, our results imply a particularly interesting 
bound in the dimension five case, T r i, > 10 9 GeV. For a reheating temperature in this range, 
the gravitino problem provides very challenging limitation on the supersymmetric parameter 
space. 

We emphasis that our results are robust in various concerns. The derived bounds do 
not depend strongly on the choice of e en d- Although our results are achieved for single¬ 
field inflation models we expect the bounds to conservatively hold for the situation in which 
inflation is ended by a second field and hence € en d is not required to become large towards 
the end of inflation. In order to match the inflaton potential at the end of inflation to the 
subsequent oscillation phase we assumed a (predominantly) quadratic potential in the latter 
phase, i.e. an equation of state ui = 0. We argued that even for a different choice of co 
our results are not changed significantly. Throughout this study we worked in the slow-roll 
approximation. Nevertheless we have discussed that our bounds are expected to hold up to an 
0(1) factor even if the slow-roll conditions are slightly violated. In contrast to the robustness 
of the bound on I4nd there, however, remains a huge dependence of the bound on TA to 
an unknown quantity, namely the type and coefficient of the operator governing the inflaton 
decay. This amounts in a difference of almost six orders of magnitude between the case of a 
dimension five and six operator. 
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Appendix 


A CMB observables beyond leading order 

In this Appendix we present the full expressions for the CMB observables in terms of the 
potential slow-roll parameters used for the computation of x 2 as described in Sec. 13.21 In all 
expressions the slow-roll parameters are evaluate at (j) = </>*. Note that there does not exist a 
slow-roll parameter beyond a 3 for the considered class of potentials (|24|) as the fifth derivative 
of the potential vanishes. To second order in the slow-roll parameters the scalar amplitude 
reads [67] 
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Using the results from [68] and m we can obtain a second order expression for the 
tensor-to-scalar ratio: 
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